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DIFFUSION  LIMIT  FOR  THE  TRANSPORT  OF 
QUASIFREE  ELECTRONS  IN  LIQUIDS 
P.  J.  Paes  Leme* 

Department  of  Chemistry,  New  York  University 
4 Washington  Place,  New  York,  NY  10003 

Abstract 

A model  of  photoelectron  emission  by  liquids  is  treated  in  terms  of  a 
transport  problem.  Its  asymptotic  behavior  is  studied  in  the  limit  when 
the  mean  free  path  becomes  small  and  the  geminetric  recombination  is 
neglected.  This  solution  is  then  shown  to  be  in  agreement  with  that 
previously  assumed  in  the  Nemec  theory. 

1.  INTRODUCTION 

1 2 

The  model  of  photoelectron  emission  by  liquids  studied  ’ involves  three 
steps:  generation  of  quasifree  electrons  in  the  bulk  of  the  liquid,  their 
transport  towards  the  liquid-vacuum  interface,  and  the  overcoming  of  the 
barrier  at  the  liquid-vacuum  interface.  The  probability  distribution  of 
finding  a quasifree  electron  generated  at  a certain  point  rQ  with  kinetic 
energy  E^,  in  a volume  element  d?  about  t with  kinetic  energy  dE  about  E, 
plays  a central  role  in  the  study  of  the  energy  distribution  curves  (EDC)  of 
photoelectrons  emitted  by  solutions.  Such  probability  distributions  can  be 
obtained  from  the  solution  of  a steady  linear  transport  equation.  We  study 

♦Present  address:  Departmento  de  Matematica,  Pontificia  Universidade  Catolica 
de  R.  J.,  R.  Marques  de  S.  Vicente  225,  22453  - Rio  de  Janeiro  - R.  J.,  Brasil 
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the  asymptotic  behavior  of  such  solution  in  the  limit  when  the  mean  free  path 

between  collisions  of  the  electron  with  the  liquid  molecules  is  small 

(diffusion  theory  limit)  and  the  Coulombic  interaction  between  the  electron 

and  its  parent  is  not  considered.  Here  we  make  use  of  known  results  from  the 

neutron  theory  of  slowing  down  problems  (Fermi-age  theory).  The  solution 

? 

we  then  obtain  agrees  in  form  with  the  solution  assumed  by  Nemec.  In  our 

case  the  law  of  energy  loss  in  the  process  can  be  derived  once  we  know  the 

form  of  the  electronic  collision  cross-section.  The  problem  treated  here,  but 

with  the  added  complication  of  geminate  recombination,  is  presently  being 

4-6 

investigated  by  adapting  different  techniques. 


We  assume: 


2.  FORMULATION 


(al)  the  electrons  are  classical  point  particles  of  mass  m,  moving  about  a 
medium  consisting  of  host  particles  of  mass  M,  with  m « M; 

(a2)  the  density  of  the  electrons  is  sufficiently  small  so  that  the  effect 
and  frequency  of  electron-electron  collisions  are  negligible; 

(a3)  the  solute  molecules  are  sufficiently  dilute  in  such  a way  that 
electron-ion  interactions  are  negligible; 

(a4)  the  liquid  is  isotropic  and  homogeneous; 

(a5)  there  is  no  absorption  (of  electrons); 

(a6)  the  coulombic  interaction  between  the  electron  and  its  parent  is  not 
considered; 


* 


Unrts*'  "tit  ai:iv/e  i unr?  (a!'— («*“'  oir  piuv:uelHCtron  e-ission  model  car 

do  fov"’.Llatec  11  C*v*T:S  of  c stead-  Vnea"  t'ar.soort  equation  for  fne  electrons 

r 

;v-“r  + vo(v)!ll(r,\)  = J v'ejv1  )f(v'-  v )N(r,v' )dv'  * S(r,v)  (2.1) 


where  N'r,v)drdv  denotes  the  number  of  electrons  in  the  element  of  volume  drdv 
about  (r,v);  v<j(v)  denotes  the  probability  per  unit  time  of  collision  of  an 
electron  whose  speed  is  v;  f(v'~  v)dv  the  probability  that  an  electron  with 
speed  v'  before  collision  achieves  a speed  in  the  element  dv  about  v;  and  S is 


a source  term.  Note  because  of  (a5) 


r 


f ( v v)dv  = 1. 


We  are  interested  in  the  solution  of  (2.1)  when  an  isotropic  source  is 
positioned  at  a point  r , which  we  can  take  to  be  the  origin.  Here  we 
assume  electrons  disappear  once  their  kinetic  energy  is  below  a certain 
specified  value.  This  value  can  be  its  thermal  energy  or  any  other  related  to 
the  interface  barrier  ( 1 iquid-vacuum)  which  is  assumed  to  be  constant. 

Let  v = vn.  The  transport  equation  (2.1)  can  be  rewritten  as  an  integral 
equation  for  the  electron  density 

n(r,v)  = j'  N(r,vn)d3  (2.2) 

To  do  this,  first  integrate  (2.1)  along  its  characteristic  curve  in  the 
following  way.  Denote  the  right-hand  side  of  (2.1)  by  Q(r,v),  since  C.  • 

•v  ■+ 

is  simply  the  derivative  taken  along  the  direction  of  w,  if  we  replace  r by 
r - Rfl  in  (2.1)  we  can  write  it  in  the  following  form 

[-  d/dR  + o(v)]N(r-Rn,vfl)  = v_1Q(r-Rn,v)  (2.3) 

after  dividing  both  sides  by  v.  If  Q(r-Rn,v)  is  supposed  known,  integration 
of  (2.3)  leads  to 


N(r,vfl)  * N(r-R  n,vn)exp  (-o(v)R  ) + 
0 R 0 


v_1  f 0 Q(r- Rfl,v)exp  (-o(v)R') 
J n 


dR ' 


(2.4) 
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after  setting  R * 0.  Because  there  are  no  electrons  coming  directly  from 
infinity  we  have 

lim  N(if-R0^,vj!)exp  ( -o(v)R0}  = 0. 

R -*OD 
0 

Therefore  we  can  write  (2.4)  as 

N(r,vn)  * v"1  Q(r-Rfl,v)exp  { -a(v)R;dR  (2.5) 

0 ...  ->  2 
If  we  integrate  (2.5)  over  all  Q,  put  r-Rf.  = r'  and  notice  that  dRdn  = dR/R 

* dr'/ir-r'|^,  we  find  for  (2.2) 

n(r,v)  * v'1  j"  dr ' ( r-r ’ |’^Q(r' ,v)exp  { -o(v) jr-r ' ] } 

or  in  terms  of  the  flux  $(r,v)  = vn(r,v) 

$(r,v)  * (4it)-^  dr ' |r-r ' | “?exp  { o(v)  |r-r ' | } • 

{ f dv'o(v' )f(v'-*  v )* ( r ' , v ' ) + 4TrS(r',v)}  (2.6) 

J0 

Equation  (2.6)  can  be  written  in  terms  of  energy  dependent  qualities  by 

2 

changing  variables,  E = mv  12.  We  get 

$(r,E)  * ( 4rr ) /dr' |r-r' | "?exp  { -o(E)  |r-r ' | ) • 

(f  dE'o(E' )f (E1-  E)  (r',E')  + 4irS(r',E)}  (2.7) 

^ ® 1/2 

where  f(v'-  v)  = (2mE 1 )iUf (E •-  E). 


3.  DIFFUSION  WITH  MODERATION  APPROXIMATION 
Assume  the  energy  after  each  collision  is,  with  equal  probability, 
anywhere  between  (l-a)E'  and  E',  for  some  fixed  value  of  ° > 0,  i.e., 
f(E'-  E)  * J 1/oE1  for  (l-a)E'  * E 4 E' 
c0  otherwise 
Then  we  can  write  (2.7)  as 
•(  r,E ) * (4t)-1  J dr ' Jr-r ' |’2exp  ( -<?(E)  |r-r'  |)  • 
rE/(l-ai) 


J, 


dE ' ( ^ ' )-1 J(E ' ) *(r' ,E ' ) + 4wS(r*,E)} 


(3.1) 
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In  terms  of  the  lethargy  u * ln(EQ/E),  where  EQ  is  so  chosen  that  the 
lethargy  would  be  positive  for  practically  all  the  electrons  which  play  a role 
in  the  theory,  for  example,  Eq  = 100  eV,  and  the  collision  density  F("r,u)  * 

Eo(E)«(r,E),  eq.  (3.1)  becomes 

F(r,u)  = (4tt)*1<j'u) /dr'  |r-r' |"2exp  { -a(u)  |r-r' | ) • 


(a-1  j du,eu,_uF(r' ,u' ) + 4nSu(r',u)) 


(3.2) 


where  ; = - 1 n( 1 - a)  and  S (r,u)  = ES(r,E)  is  the  number  of  electrons 

liberated  per  unit  volume  and  unit  time  and  per  unit  lethargy. 

Equation  (3.2)  can  be  transformed  by  Taylor's  expansion,  into  a 
differential  equation  if  we  assume  S and  F do  not  vary  much  within  a mean  free 


path  l/o.  Expand  F in  powers  of  (u'-u),  then 

F(r',u')  = F(r',u)  - (u-u1 ) ( 3/3u)F(r ' ,u)  + higher  order  terms 

Because 

a'1  /°  du'eu  ’u  = 1 
u-c 
and 

a*1  f du'(u-u')eu  u * a *{-  ce’C-  e‘C+  1) 
u-c  i 

= 1 + a 1(l-a)ln(l-cx)  = £ 

(£  * average  increase  in  u),  from  (3.3)  we  obtain 


(3.3) 


-1  fU  du'eu"uF(r'  ,u' ) = F(r',u)  - ( a/3u)  CF(r  • u) 


(3.4) 


if  we  neglect  the  higher  order  terms.  The  variation  of  F with  respect  to  u is 
caused  by  its  variation  with  r and  it  will  be  small  if  the  latter  is  also 
small.  This  is  a condition  for  the  validity  of  this  development. 

Expand  F(r',u)  in  powers  of  (r-r* ), 

F(r',u)  » F(r,u)  ♦ \ (r-r1 ) i ( 3/3ri )F(r,u) 


+ 2-il  ! 


i=i  j=i 


(r-r,)i(r-r’)j(32/3ri3rj)F(-jU) 


♦ higher  order  terms 


(3.5) 
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where  the  subscript  i stands  for  the  i-th  component.  Because  of  the  factor 
exp  { -o(u) 1 r-f ' | } in  (3.2),  only  the  contribution  for  small  values  of  |~r-"r  • j , 
up  to  the  order  of  a few  mean  free  paths  1/  a(u),  will  be  significant. 
Therefore  we  replace  F(r’,u)  in  the  first  term  of  (3.4)  by  the  first  three 
terms  of  the  series  (3.5),  Su(r',u)  and  (3/3u)£F(r* ,u)  by  their  values  at  r' 
= r.  Since 

(4tt)_1j  J dr'  J r-r 1 |"2exp  { -cr|r-?'  J } = 1 , 

(4tt)'*7  Jdr*  | r-r 1 j ~2( r-r ’ ) ^exp  { -a|r-?'  J } = 0 (i  = 1,2,3) 
and 

(4^)  J"  dr ' | r-r  * |”2(r-r' ) i (r-r*  J^.exp  { -a|r-r‘  |>  = (2/3a2)5^ 

we  get  from  (3.2) 

(3/3uUF(r,u)  = (l/3c2(u))  AfF(r,u)  + Su(?,u)  (3.6) 

4.  COMPARISON  WITH  NEMEC'S  APPROXIMATION 
Introduce  a new  variable  t,  referred  to  as  the  age  of  the  electron, 

r=  fU  du/35o2(u)  = f E°  dE/3cEo2(E)  (4.1) 

\ E 

The  diffusion  equation  (3.6)  then  becomes 

(3/3t)(5F)  = Ap(5F)  + ST(r,T)  (4.2) 

where  S T ( r , t ) = 35a  (u)S  (r,u)  is  the  number  of  electrons  produced  by  the 
source  per  unit  time  and  unit  t interval. 

For  a point  source  of  monoenergetic  electrons  and  unit  strength,  equation 
(4.2)  gives  the  distribution 

q = CF  = [4ir(x-Tg)  ]"2/^2exp  { -r2/4(i-tg)}  (4.3) 

where  is  the  age  of  the  electrons  emitted  by  the  source.  If  the  source 
is  not  monochromatic,  the  solution  of  (4.2)  is  obtained  by  integrating  (4.3) 
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against  the  distribution  function  of  the  energies  or  ages  of  the  generated 
electrons. 

2 

The  solution  (4.3)  here  obtained  agrees  with  Nemec's  if  we  identify  the 

1/2 

quantity  he  calls  R = (i.X/6)  , where  I is  the  total  length  of  the  random 

walk  and  X is  its  mean  free  path,  with  the  square  root  of  t-t  , i.e.,  if 

R2  = t - t = f 9 dE/3SEa2(E)  (4.4) 

9 

If  we  assume  a(E)  = o where  cq  is  a constant,  then  from  (4.4)  we  get 


R = (3£a2f1/2;in(Eg/E)]1/2 
If  o(E)  = °0E'n/2,  n * 0,  then 
R = (35a2) "1/2(Eg  - En)1/2 

The  laws  of  energy  loss  considered  by  Nemec  can  be  recovered  from  the  above 
for  n = 3/2,  2 and  1 
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